Quantum simulation of Anderson and Kondo lattices with superconducting qubits 
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We introduce a mapping between a variety of superconducting circuits and a family of Hamil- 
tonians describing localized magnetic impurities interacting with conduction bands. This includes 
the Anderson model, the single impurity one- and two-channel Kondo problem, as well as the ID 
Kondo lattice. We compare the requirements for performing quantum simulations using the pro- 
posed circuits to those of universal quantum computation with superconducting qubits, singling out 
the specific challenges that will have to be addressed. 



I. INTRODUCTION 

Quantum simulation^ consists on tuning the dynamics 
of a flexible quantum mechanical system to simulate the 
properties of another physical system or of a quantum 
mechanical model whose solution is unknown. While a 
universal quantum computer can efficiently simulate the 
dynamics of any quantum systerrPEl, designing a quan- 
tum simulator for a specific model should simplify the 
experimental requirements. In particular, quantum sim- 
ulators require neither high fidelity thresholds nor error 
correction, and thus fewer qubits are needed to get inter- 
esting results which are out of reach for classical numeri- 
cal computations. The successful simulation of the Bose- 
Hubbard model using cold atoms in optical lattices^ is 
a paradigmatic example. 

In the last years, there have been tremendous ex- 
perimental achievements in the context of supercon- 
ducting qubits, improving both qu antum control and 
coherence time d 6 ' 7 ' WU | H | ^ | i W5 | 

While the aim is 

the implementation of a scalable universal quantum 



computer 



141151161 



this technology can be used in the quan- 



tum simulatio n of ma ny-body physics with Josephson- 
Junction arrays^ 18 * 19 ! 

With this motivation, we first design a circuit to sim- 
ulate the Anderson model for an impurity ion coupled to 
conduction electrons^!. This is done through an exact 
mapping from the effective models describing the super- 
conducting circuits to the fcrmionic Hamiltonian for the 
impurity problem. In a certain parameter regime of our 
circuit, this leads to the Kondo Hamiltonian modeling 
the non-trivial physics of itinerant electrons that inter- 
act with a localized impurity having non-zero magnetic 
momeniPS. This is one of the central models for strongly 
correlated electrons with implications going beyond Con- 
densed Matter Physics. Furthermore, our constructions 
are versatile enough to simulate other relevant Hamil- 
tonians. With an appropriate choice of the geometry, 
these include the Kondo lattice^ and the multichannel 
Kondo modeP^l The latter exhibits the simplest example 
of non- Fermi liquid behavior, which is a strong deviation 
from the standard model of condensed matter systems 
introduced by Landau. 

We will also discuss the experimental challenges for im- 



plementing such circuits. In particular, while coherence 
and interaction times of current superconducting charge 
qubits seem good enough to simulate the long-time dy- 
namical and static properties of these lattices, further 
experimental work is needed. Finally, we propose a set 
of measurements to gather information about transport, 
correlation and energy spectra and comment possible im- 
plications of the present work. 



II. CIRCUIT DESIGN 
A. Elementary components 

Let us start by describing the proposed architecture. 
We consider two or more lines of low-capacitance su- 
perconducting islands coupled by Josephson junctions 
[Fig. [Tja.] . While not shown here, each island is itself 
part of a circuit like that of a charge qubit, with an ex- 
ternal lead which capacitively induces an offset potential 
on the island, and possibly a Josephson coupling to a 
superconducting reservoir of Cooper pairs^. 

Our quantum simulation protocols are based on two 
building blocks. First, as we will show below, each island 
will have a very low capacitance and be constrained to 
have at most one Cooper pair, which can be treated as 
an impenetrable boson. Second, selected pairs of islands 
will be coupled capacitively [Fig. TJ)]. The coupling and 
voltages of these islands will be adjusted to create a suit- 
able energy landscape that favors having a single excess 
Cooper pair in either the upper or lower island, states 
which we associate with an effective pseudo spin. The 
appropriate choice of energies has the following form 



E(n ^n ol ) = Un ^n Ql - e(n 0T + n ;). 



(1) 



Here, nof and no^ represent the excess of Cooper pairs 
on the upper and lower island in the impurity [Fig.[lk-b]. 
The constants U and e denote, respectively, the repulsive 
interaction energy coming from the capacitive coupling 
between islands and an energy offset of these islands with 
respect to the rest of the circuit. 

As mentioned before, when we embed such an element 
in a circuit, the energy offset e enforces the impurity to 
host at most a single excess Cooper pair. In Fig. [TJ: we 
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FIG. 1: (a) Quantum circuit simulating the single Anderson 
or Kondo impurity model. The big square boxes are super- 
conducting islands, connected by crossed squares which are 
Josephson junctions. A selected pair of islands is capacitively 
coupled and represents a magnetic impurity. The Josephson 
energy for tunneling between junctions are J for the main 
lattice and t around the impurity, (b) Two occupation states 
for the impurity representing a pseudo-spin state, (c) Energy 
levels of the impurity and a neighboring site, not, n oi an d N r 
represent the number of pairs of the two islands in the impu- 
rity and of the rest of the circuit. By external voltages, — e/2, 
and the capacitive coupling, U ~ 1/C, it becomes favorable 
to have a single pair in the impurity. 



have depicted the lowest energy levels of a circuit con- 
taining a single impurity, with occupation numbers not 
and no i, connected to an array of ordinary islands, with 
total population N r . A possible configuration is shown in 
Fig.[lji. If wc neglect the kinetic energy of hopping pairs, 
we find that the states where the impurity is occupied by 
a single pair, rt | = 1 or noj, = 1, have lower energy than 
states where this pair is in the rest of the circuit. Multi- 
ply occupied impurities, n T + n o| > 1, have a even larger 
energy, U, so that we can neglect these states. 



B. Quantum circuit model 



physics of a Josephson junction array such as the one in 
Fig. Ei 

H = £ \{n ia - n^Ui^n^ - n$) 

+ X! JijCos((f>i a ~ cf) ja ) + E(n Qh n ol ). (2) 

Here, n irT and ($>i a are conjugate variables denoting the 
excess number of Cooper pairs on an island and the as- 
sociated superconducting phase, respectively. is the 
Josephson energy for a junction that connects two is- 
lands. Finally, there is a short range interaction term, 
U = 4e 2 C~ 1 , which is the inverse of the capacitance ma- 
trix. Similar to J, C is a sparse matrix with nonzero ele- 
ments on the diagonal, denoting the intrinsic capacitance 
of each island, Ca ^ 0, and also between islands con- 
nected by junctions or capacitors. We choose all capaci- 
tances to be small, 4e 2 /C 3> J, and impose Ca 3> Q^j 
so that we can neglect any off-site interaction, U^j ~ 0. 

The residual charges, nf — n- + n\ D , are a combi- 
nation of a possible systematic shift, nf^ , and an effective 
( DC ^ — n(s)-\/-{9) io c proportional to the potential 



term n)""' = C} 9> V} 3 ' /2e, 

r{9) 



applied to each island, V i , and the capacitance through 



which it is applied, C\ y> [See Ref. [24] From this it follows 
that the residual charges can be set to zero by appropri- 
ate tuning the offset voltages acting on eac h junction. As 
in current single and two-qubit experiment d 7 * 8 * 11 * 12 * 14 ^ we 
will work in the so called sweet spot of energetic degen- 
eracy, given by — 1/2. Due to the low capacitance, 
U 3> J, the state of each island is then constrained to 
having or 1 Cooper pair, or any linear superposition of 
these states. 

Contrary to the cavity-QED and superconducting 
qubit experiments^! we would like the total number of 
Cooper pairs to be a well defined quantum number. This 
means we do not want a residual coupling between the 
states with and 1 pairs, which would correspond to a 
term Ja cos(0;) in the previous model. This is effectively 
achieved by switching off the coupling between the super- 
conducting islands and the Cooper pair reservoir during 
the experiment. 

Finally, it is interesting to note that the total num- 
ber of pairs in the ground state can be tuned with very 
small voltages that take islands away from degeneracy 
and which play the role of a local chemical potential 



t'i, 



1 1 " 



(3) 



In Sect. |II D| through |II G] we will combine the previ- 
ous elements in circuits that simulate specific Hamiltoni- 
ans. In order to perform this task we must first introduce 
the mathematical description of such circuits and the ap- 
proximations under which they will oper ate. We start 
considering the standard Hamiltonia n 118 ! 1 ^ describing the 



C. Hard-core boson model 

Summarizing the previous approximations, we have an 
array of superconducting islands coupled both capaci- 
tively and via Josephson tunneling. At the same time 
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we have ensured that each superconducting island has 
o 1 excess Cooper pair, all other states being energeti- 
cally unfavorable. In other words, we have imposed the 
hard-core bosons conditions on the Cooper pairs: they 
become like impenetrable bosonic particles such that no 
two pairs with the same effective "spin"^ can coexist on 
the same site. 

In this limit we can describe any configuration of the 
circuit using the occupation numbers {n^}. Due to the 
hard-core condition, the dynamics is well approximated 
by the projection of the circuit Hamiltonian Q onto the 
relevant energy subspace, ni a £ {0, 1}. To perform this 
projection we will need to express the phase terms in 
basis of occupation numbers 

y~] Jij cos{^ ia - cj) ja ) = (4) 

ija 

/A^ii ^ \riia + l,nj <T -l){n i<T ,nj <T \+'H..c. 

Note that this equation has a very simple interpretation: 
Cooper pairs are transferred one by one between any two 
islands connected by the hopping matrix Jij. 

The projection is implicitly performed by means of a 
more concise mathematical representation that includes 
an exclusion principle for particles on the same site. We 
introduce hard-core bosonic operators, a\ a and a^, which 
respectively create and annihilate a Cooper pair on a 
given site. They act on the basis elements as follows 

Oio-Ko- = 1) = K<t = 0), a ia \0)=0, 



A 



a ia \u ) = \l), at|l) = 0, 
have the usual bosonic commutation relations 



(5) 



and enforce the hard-core condition simply because = 
0. Using these operators we can express the number of 
particles n.; CT = a\ a a,i a and the projected phase operator, 

2cos((/) J : cr - (f>ja) ^ a\ a a jc , + a} ]a a iai (6) 
which adopts the form of a hopping term. 

D. Single impurity Anderson model 



Out of all terms in the hopping matrix J^, we have left 
only the Josephson energies between sites in the conduc- 
tion array, Ju+i = J for i > 0, and the coupling be- 
tween the impurity site and this band, Jqi = Jio = t 
[See Fig. [lji] . Finally, as mentioned in Sect. 



II C we have 



included an effective chemical potential which depends 
on deviations from the degeneracy point of the local po- 
tentials acting on each island. 

We can map this Hamiltonian to a fcrmionic 
model using a nothe r standard tool, the Jordan- Wigner 
transformatiorPSEH 



a 4T = ^(-IjEi^t, 
an = c u (-l)^<^ +J \ 



(8) 



c ia C io ~ a icr a icr- 



The new operators satisfy the usual fermionic anti- 
commutation relations {ci a ,Cj a i} — 5ijS a a' and N a — 
'^2i n i(T- With this transformation, our effective Hamilto- 
nian (JsJ) converts identically into a fermionic model 

H A = -J^(ct +lCT c iCT +H.c)-t^(cJ ff c lff +H.c.) 

i>l a 

+ E(n oh n ol ) +^2/J, ir7 n i(r , (9) 

itj 

which is known as the Anderson impurity modeP^. We 
have thus demonstrated that the Anderson Hamiltonian 
accurately describes the dynamics of a circuit that we 
have designed, under reasonable approximations. Fol- 
lowing the ideas sketched in the introduction, we now 
argue that the quantum circuit can itself be used to sim- 
ulate the dynamical and static properties of the Anderson 
model. For instance, if we are able to create and cool this 
circuit to its ground state, we will create a state which is 
related to the ground state of the Anderson model by a 
unitary transformation. The properties of the equivalent 
fcrmionic system can be recovered from measurements of 
the Cooper pairs and the Jordan- Wigner relations 

In the following subsections we will show how, by 
choosing other limiting cases and designing slightly dif- 
ferent circuits, one can simulate not only the Ander- 
son model, but other more complicated and interesting 
Hamiltonians. 



We will apply all the approximations and techniques 
introduced before to study a particular circuit, which 
consists of two impurity islands coupled to two ID 
Josephson junction arrays [Fig. [Ik]. The low energy dy- 
namics of this circuit is described by the bosonic Hamil- 
tonian 

H A = -J^(aJ +lCT a iCT +H.c)-t^(aJ (r ai CT +H.c.) 

i>l a 

+ E(n oh n ol ) +^2^i ia n ia . (7) 

ia 



E. Kondo model 

A special limit of the Anderson model is obtained when 
we impose a low capacitance on the impurity, U — > oo, 
and make the energy displacement large compared to the 
coupling between the impurity and its neighbors, t <C 
e. Using standard second order perturbation theory, the 
result is the single impurity Kondo HamiltoniarP^l 

H K = ~JJ2 c i<r°i° + Jk $o ■ si- (10) 
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FIG. 2: (a) Circuit for a two-channel single impurity Kondo 
system. Each half of the circuit is associated to a different 
channel. An extra junction with negative sign, t' ~ —Jk, is a 
quantum interference element that prevents hopping between 
channels, (b) Circuit for a ID Kondo lattice. 



Here So and si are the pseudo-spins of the impurity and 
of the first lattice site 



So 



a.a 1 



' C L C 0> 



tj' ) and s i = ^(Tqo-'C^cio-'^II) 



c7,<t' 



respectively, and = t 2 /e is the anti- ferromagnetic cou- 
pling between the Kondo impurity and the free fermions. 
We want to stress that by tuning the parameters J, t and 
e, we can realize a crossover from a low energy regime be- 
low the Kondo temperature Tk to a high energy limit. 



F. Two-channel Kondo model 

The previous setup and treatments can be reused to 
simulate a two-channel Kondo systenPH. As shown in 
Fig. [2^i, the impurity is in the middle of a one-dimensional 
lattice, coupled to two different conduction bands. Fol- 
lowing similar steps, the new effective model becomes 

H K = - J J3 c Ura C 3<r<* + J K 53 ^° ' * lc " ( 12 ) 
(i,j)(7,a a — 1,2 

where a S {1,2} denotes the two fermionic channels of 
our problem. Note that in this Hamiltonian we have 



removed the terms that induce hopping between different 
channels, JKc\ aa iC\ aa . These and higher order couplings 
are canceled by two cleverly placed small junctions whose 
tunneling has been magnetically calibrated to the value 
t'~-J K . 



G. Kondo lattice 

A slightly less trivial circuit is required to simulate 
a Kondo lattice that has multiple localized impurities 
interacting with fermions. There are two different ways 
of doing it, leading to slightly different physics in each 
case. A simple concatenation of the circuit in Fig. [2Jl will 
produces a set of impurities coupled by a number of sites 
with free fermions, equivalent to the model introduced by 
Paredes^ZI for cold atoms. However, the Kondo lattice as 
known in the literature^ considers multiple impurities 
coupled to a common conduction band, which can be 
achieved using the circuit in Fig. [2Jd. Following Sect. |II C| 
we can write down an effective bosonic Hamiltonian 



H = -J^a| a a j(T -t^(& t feff a. ifcff + H.c) 

+ ££(4?-4Y)+E^ ( - 

k iu 



(13) 



where we have introduced different hard-core bosonic op- 
erators for pairs in the conduction band, denoted by 
ai a ,a\ a , and — a^dia, and for pairs in the impu- 
rities, given by &i<r,&L, and njjfj = b\ a b ka . 

This circuit does not have a direct translation into a 
fermionic Anderson model because for such a quasi-2D 
structure there is no Jordan- Wigner transformation that 
preserves the shape of the bosonic Hamiltonian. We will 
however perform an incomplete Jordan- Wigner transfor- 
mation which acts only on the conduction arrays, a irT1 
leaving the bka operators untouched, 



an = c a (-i)^-^'+< ) ) 



(14) 



where — ^ n \a ■ With this unitary map, and work- 
ing out the second order perturbation theory, one arrives 
to the standard Kondo lattice model 



H K L = -J E C i<? C 3 cr+ 53 ~Y Sk '°aa'c[ a C lc r>. (15) 



Here, Sk denotes the spin of the fc-th impurity which is 
coupled to the site i% and is defined using bosonic oper- 
ators Si k = J2 a .a> \^aa'b\ k0 b lk a- ■ This general notation 
allows for having less impurities than conduction elec- 
tron, as well as them being placed on arbitrary positions. 
The coupling is once more anti- ferromagnetic Jk ~ t 2 /e 
and just as tunable as in the case of a single impurity. 
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III. ENERGY SCALES AND DECOHERENCE 



Let us consider what are the experimental challenges 
for implementing these mappings. We begin with the 
energy scales of our problem. On the one hand we have 
free fermions, which hop with a tunneling amplitude J 
and a Fermi energy ep ~ 2 J. In current experiments with 
charge qubits J can be roughly of order of a few gigahertz. 
On the other hand we have the Kondo coupling Jk — 
t 2 /e, which will be in general smaller. For this coupling 
to be relevant, the temperature at which the experiment 
is performed should lay below the Kondo Temperature, 
^bTk = £f exp[l/ Jjcpikp)], where p(kf) is the density of 
states at the Fermi energy. Since in our case the density 
of states at half filling (fi = 0) is roughly p(kf) = n/J, 
using t ~ J, we have UbTk ~ 2Jexp(— e/Jir). This 
means that, conservatively, for e/J — 8 — 10, the Kondo 
temperature Tk ~ 40 — 75mK would be larger 20 mK, 
the refrigerator temperature of current experiments^. 

Another important challenge is decoherence, which in 
our setup arises mainly from charge fluctuations. As is 
the case of scalable quantum computation, the charac- 
terization of decoherence in multi-qubit setups is still 
an open problem that deserves further investigation. As 
a guide, we consider the time-s cale f or the decoherence 
of a single qubit is T2 ~ which is achieved 

for dispersive readouts. This time-scale gives a fre- 
quency of 2 MHz which is well below J/h ~ 10GHz and 
two orders of magnitude away of the lowest frequency 
Jic/h = t 2 / e ~ 100 MHz that we find here. Furthermore, 
associated to this decoherence rate we can establish an 
effective temperature of 0.1 mK which is also well below 
the Kondo temperature. 

When compared to other systems, such as cold atoms 
in optical lattices^, the superconducting circuit approach 
looks as enjoying real advantages. First, the supercon- 
ducting qubits are not constrained to periodic or quasi- 
periodic structures. Second and most important, the in- 
troduced setup has the potential quality of a good quan- 
tum simulator, which is measured by the number of times 
a particle can hop before the wave function is affected 
by decoherence. For cold atoms in the lowest band of 
an optical lattice^, coherence times of 500 ms are to be 
compared with a hopping rate of 1kHz between lattice 
sites, giving a ratio of 500 or better. A similar ratio with 
superconducting circuits, where the hopping is around 10 
GHz, means the decoherence time can be 0.05 us, an order 
of magnitude faster than current experiment er ^ 28 !. On the 
other hand, superconducting circuits have the disadvan- 
tage of being one-time experiments because the geometry 
of the couplings cannot be modified in real-time. In ad- 
dition, it may be complicated to calibrate all islands to 
reach the sweet spot. However, if we realize that the 
couplings can be initially switched to zero, this task is as 
difficult as tuning the individual superconducting qubits 
in a scalable quantum computer. 



IV. MEASUREMENT 

While we have demonstrated that certain quantum cir- 
cuits can be used to simulate interesting fermionic mod- 
els, the equivalence between circuit and the Hamiltonian 
involves the use of Jordan- Wigner transformations pj. 
Therefore, not all measurements on the quantum circuit 
will give directly properties of the respective fermionic 
system. One can envision three kinds of measurements in 
these systems which do not have this problem. The first 
one are global transport measurements that mimic past 
experiments with Josephson junction arrays^. Since 
they are based on total particle numbers, they are the 
same for fermionic and bosonic models. 

Next we have spectroscopic measurements and mea- 
surement of energy gaps. The coupling between the im- 
purity and the conduction band leads to a binding en- 
ergy of order A = ksTfc given above. This energy can 
be measured by trying to polarize the impurity with an 
oscillating voltage applied to the island. We expect this 
to have an effect either when the potential is of order A 
or when it oscillates with a frequency of order A/h. 

The third kind of experiments corresponds to mea- 
suring the individual qubits and extracting information 
about certain correlation functions. In the case of Fig. [pi, 
the formation of a singlet with the conduction band elec- 
trons gives rise to a strong spin correlation in the vicin- 
ity of the impurity. The correlation length is of order 
£ = exp[l/ Jxpikp)} sites and can be seen in the con- 
nected correlator as a function of the site m 

C m = (SozS mz ) - (S 0z )(s mz ) ~ exp(-m/£). (16) 

This is known as the Kondo cloud and has not been di- 
rectly detected by condensed matter setups. In our case, 
C m can be computed from the statistics of the population 
difference between islands, Sj z = (njt — flit ) /2, which is 
best measured using dispersive probe^SI22l 



V. CONCLUSIONS 

Summing up, in this paper we have presented a map- 
ping between some superconducting circuits to a family 
of relevant Hamiltonians describing magnetic impurities 
in conduction bands. While there are still implementa- 
tion issues to be solved, they seem less demanding than 
those of a fully scalable quantum computer implementa- 
tion with the same technology. We thus expect that our 
work will motivate further developments, both on the ex- 
perimental and on the theoretical side, searching better 
multi-qubit coherence times and other mappings which 
are based on flux or hybrid qubits. 
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